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Abstract: The maximum mass of a neutron star (NS) is poorly defined. Theoretical attempts to define this 
mass have thus far been unsuccessful. Observational results currently provide the only means of narrowing this 
mass range down. Eclipsing X-ray binary (XRB) pulsar systems are the only interacting binaries in which the 
mass of the NS may be measured directly. Only 10 such systems are known to exist, 6 of which have yielded 
NS masses in the range 1.06 - 1.86 M .We present the first orbital solutions of two further eclipsing systems, 
OAO 1657-415 and EXO 1722-363, whose donor stars have only recently been identified. Using observations 
obtained using the VLT/ISAAC NIR spectrograph, our initial work was concerned with providing an accurate 
spectral classification of the two counterpart stars, leading to a consistent explanation of the mechanism for spin 
period evolution of OAO 1657-415. Calculating radial velocities allowed orbital solutions for both systems to 
be computed. These are the first accurate determinations of the NS and counterpart masses in XRB pulsar 
systems to be made employing NIR spectroscopy. 

1 Introduction 

Despite extensive and ongoing theoretical work on the NS equation of state (EOS), the precise nature 
of the fundamental physical properties of NS matter is still poorly defined. Observational work can 
assist in reducing the number of contending theories by eliminating those that place unrealistic con- 
straints on the mass range of observed NS. NS masses may only be determined from binary systems, 
within this paper we consider a specific class of these objects, those containing an eclipsing X-ray 
binary pulsar. At present only 10 such systems are known within our Galaxy. Prior to this work 
6 of the NS in these systems have had mass determinations, with the donor star in each observable 
optically. Within this paper we discuss the first mass estimates found for NS employing near-infrared 
(NIR) spectroscopy. We have studied two eclipsing X-ray pulsar systems containing a High Mass 
donor, OAO 1657-415 and EXO 1722-363. Initially an accurate spectral classification for each of the 
donor stars was conducted utilising observations made using the VLT/ISAAC NIR spectrograph and 
current NIR spectral atlases. Using multi-epoch NIR spectra of each system we were able to deter- 
mine the radial velocity of the donor star in each of the two systems thus enabling the construction of 
an orbital solution. This solution was then employed in calculating the mass estimate of each NS and 
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the corresponding high mass donor, placing constraints upon the system inclination and separation of 
the binary system. 

2 Spectral classification 

2.1 Spectral classification of EXO 1722-363 

EXO 1722-363 (alternatively designated IGR J 17252-36 16) was discovered in 1984 by EXOSAT 
Galactic plane observations (Warwick et al. 1988). XMM-Newton observations narrowed down 
the source position location to 4". This allowed the identification of an IR counterpart 2MASS 
J17251 139-3616575 (with magnitude J = 14.2, H = 11.8 and K s = 10.7 (Zurita Heras et al. 2006)). 
Examining Fig.l we can see that all of absorption lines in this spectrum are narrow, indicative of the 
object being a supergiant. EXO 1722-363 shows the singlet He I 2.058 fim line in emission, this line 
being highly sensitive to wind and temperature properties. The N ill 2. 1 15 /im emission line is a com- 
mon feature in B0-B1 supergiants. The absence of strong Br7 2.1655 /im emission features implies 
that EXO 1722-363 does not exhibit a strong stellar wind. From a qualitative comparison of spectra 
from Hanson et al. 2005, we identify EXO 1722-363 as being of spectral type B0-B1 la (Mason 
et al. 2009). By comparison with evolutionary rotational massive star models (Meynet & Maeder, 
2000) we find an initial progenitor mass for EXO 1722-363 in the range 30M o - 4OM . Following 
the method for determining spectroscopic distance as detailed in Bibby et al. 2008, we determined a 
distance for EXO 1722-363 of 8.0t|o kpc which is comparable within errors to the distance deduced 
in Thompson et al, 2007. Comparing our calculated distance with model fluxes derived from spectral 
fits to EXO 1722-363 (Corbet et al, 2005), we found that EXO 1722-363 has an intrinsic X-ray flux 
variability (in the range 2-60 keV) such that F min = 0.78 x 10~ 10 erg cm~ 2 s _1 and F max = 12.2 x 
10~ 10 erg cm~ 2 s -1 . We derive X-ray luminosities for EXO 1722-363 such that L Xmin = 3.4 x 10 35 
erg s _1 and L Xmax = 1.6 x 10 37 erg s -1 . We find this luminosity range entirely consistent with EXO 
1722-363 being the donor within an SGXRB system. 

2.2 Spectral classification of OAO 1657-415 

OAO 1657-415 was first detected over 30 years ago by the Copernicus X-ray satellite. From an 
examination of the orbital parameters of the X-ray pulsar it was determined that OAO 1657-415 was 
a high-mass system, indicating the mass of the donor lay between 14-18 M with a corresponding 
radius range of 25-32R . Determination of these stellar parameters led to a suggested classification 
of BO-6 I (Chakrabarty et al, 1993).The correct identification of the donor in this system required 
the precise location of OAO 1657-415 to be accurately made. This was achieved by the Chandra 
X-Ray Observatory narrowing the X-ray location error radius down to 0.5". Optical imaging of this 
position did not detect any donor candidates down to a magnitude of V>23. Near infrared imaging 
was employed to overcome significant levels of interstellar reddening, resulting in the identification 
of a donor located within the Chandra error radius. A corresponding IR counterpart was located in 
the 2MASS catalogue, 2MASS J17004888-4139214 (with magnitudes J = 14.1, H = 11.7 and K s = 
10.4) with A y = 20.4 ± 1.3, located at a distance of 6.4 ± 1.5 kpc (Chakrabarty et al, 2002). NIR 
K s band spectroscopy of the donor obtained in 2008 (Mason et al, 2009) led to a re-evaluation of 
the spectral classification. Close examination revealed that OAO 1657-415 shared a similar spectral 
morphology with that of Ofpe/WNL stars. These are stars in transition between the OB main sequence 
and hydrogen depleted Wolf-Rayet stars, whose evolution follows from a wide range of progenitor 
masses. The spectrum of the mass donor in OA01657-415 is presented in Fig. 1, and is dominated by 
He I 2.058 /im and Br7 emission, the former stronger than the latter. We find a poor correspondence 
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Figure 1: Left Figure : Comparison of EXO 1722-363 and template 09-B3 la spectra from Hanson 
et al., 2005. Right Figure: Topmost shows a spectrum of OA01657-415 compared with a Ofpe/WNL 
spectrum of IRS16NW from Martins et al., 2007. 

with the spectra of BO-6 supergiants (Hanson et al, 1996, Hanson et al, 2005) - as suggested for 
the mass donor by Chakrabarty et al, 2002 on the basis of a combination of photometric and X-ray 
data. However, comparison with the spectra of massive transitional objects presented by Morris et 
al, 1996 is more encouraging. In particular OAO 1657-415 shows pronounced similarities to the hot 
Ofpe/WNL stars. Consequently we may not a priori determine a unique distance to OAO 1657-415 
on the basis of this classification. We thus find inevitably unconstructive limits of 4.4 kpc < d < 
12 kpc. In turn this results in 1.5 x 10 36 erg -1 < Lx < 10 37 erg s _1 , also entirely consistent 
with observed luminosities of SGXRBs. Adopting the distance derived by Audley et al, 2006 leads 
to log(L/L ) ~ 5.7. For such a luminosity, comparison to the evolutionary tracks for massive stars 
(Meynet & Maeder, 2000) imply an initial mass of ~40 M . 

3 OAO 1657-415 : A mechanism for spin-period evolution 

We now turn to the implications of the Ofpe/WNL classification for the X-ray properties of OAO 
1657-415. The anomalous position of OAO 1657-415 within the Corbet diagram (Figure 2) (Corbet 
et al, 1986), is then naturally explained in terms of the properties of its stellar wind. Compared to 
normal OB supergiants (Crowther et al, 2006), Ofpe/WNL stars typically demonstrate systematically 
lower wind velocities and higher mass loss rates (Martins et al. 2007). This combination of wind 
properties permits a higher accretion rate and hence transfer of angular momentum to the NS, in turn 



leading to a smaller (instantaneous) equilibrium spin period with respect to normal OB supergiants 
(P spin oc M-^v 1 ^ 7 from Eqn. 12 of Waters et al, 1989, where P spin , M and Voo are the spin period 
of the NS and the mass loss rate and terminal velocity of the mass donor wind respectively). 
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Figure 2: Corbet diagram marking position of OAO 1657-415 and other HMXBs. OAO 1657-415 
is marked by an X, EXO 1722-363 by a filled diamond. SGXRB Roche-Lobe Overflow systems 
(Squares), Be/X binaries (Triangles), SGXRB Wind-fed systems (Diamonds) and anomalous systems 
(+) 



4 Orbital solution for EXO 1722-363 

The orbital solution we have calculated was obtained from archival ESO VLT data. Using a small 
subset of the available archive data, (11 spectra taken at different epochs spanning a wide range of 
orbital phase, from a set of 104 in total) we were able to measure radial velocities and construct the 
orbital solution shown in Fig. 3 (left). These spectra were centered on 2.1 fim, having an integration 
time of 700s, and were obtained using the SW MRes mode with a 0.6" slit. This resulted in high S/N 
spectra at a resolution R p=* 4200. The resulting NS mass that we have determined from our orbital 
solution for EXO 1722-363 is consistent with the canonical mass of 1.4 M measured in most other 
eclipsing HMXBs, except for that in Vela X-l, (Quaintrell et al., 2003). The NS mass range we have 
determined stems from a lower and upper limit obtained using the following constraints - Lower : the 
system is viewed edge on (i.e. i = 90°), Upper : the donor star fills its Roche lobe. Utilising this orbital 
solution we find a NS mass range of 1.5 - 1.6 M (Mason et al. 2010). In a similar way the measured 
mass and radius of the supergiant donor, M ~ 13 — 15 M and R ~ 25 — 28 R is determined, and 
this lends support to the B0-1 la spectral classification that we previously found (Mason et al., 2009). 



5 Orbital solution for OAO 1657-415 

As the mass donor in OAO 1657-415 is faint (H ~ 1 1.7) we employed the NIR spectrometer ISAAC 
on the VLT to obtain high resolution (R ~ 3000) and S/N spectra in the H band. Observations were 
conducted between 2008 May 13th and 2008 September 25th in the SW MRes mode with a 0.8" slit. 
Cross-correlation was performed using the standard IRAF routine fxcor. 12 high quality spectra were 
obtained that covered a wide range of orbital phase, sufficient to determine a dynamical mass solution 
for OAO 1657-415 (Fig. 3). Utilising this orbital solution we find a NS mass range of ps 1.4 - 1.7 M 
with a corresponding mass range for the counterpart star of ps 14 - 17 M . For a more precise mass 
determination please refer to Mason et al, 201 1. 




Figure 3: Left : Radial velocity data for the donor in EXO 1722-363. Right : Radial velocity data 
for the donor in OAO 1657-415. In both cases the solid line is the best fitting sinusoid with three free 
parameters, the dashed line is that with a fixed zero phase in line with the published ephemeris. In the 
case of EXO 1722-363 the orbital phase is based upon the ephemeris of Thompson et al, 2007. For 
OAO 1657-415 the orbital phase is based upon the ephemeris of Bildsten et al, 1997. 
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